Aging and domain growth in the two— dimensional 
Ising spin glass model 



H. Rieger, B. Steckemetz and M. Schreckenberg 
Institut fur Theoretische Physik 
Universitdt zu Koln 
50937 Koln, Germany 
(April 18, 1994) 

Abstract 

Interrupted aging in the two-dimensional Ising spin glass model with Gaus- 
sian couplings is established and investigated via extensive Monte-Carlo sim- 
ulations. The spin autocorrelation function scales with t/r(t w ), where t w 
is the waiting time and r is equal to t w for waiting times smaller than the 
equilibration time r eq . The spatial correlations scale with r/£(t w ), where the 
correlation length £ gives information about the averaged domain size in the 
system. Our results are better compatible with an algebraic growth law for 

£(t w ), although it can also nicely be fitted to (logi^) 1 /^ with ijj ~ 0.63. 
75.10N, 75.50L, 75.40G. 
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Aging is a characteristic feature of the non-equilibrium dynamics of spin glasses [|TJ. 
After the system has been quenched into a typical non-equilibrium state (for instance after 
decreasing the temperature rapidly) the huge relaxation times at low temperatures lead to 
the experimental observation of history dependent phenomena on laboratory time scales 
0. This scenario, also reported for various other strongly disordered materials is not 
connected to a finite-temperature spin glass transition It appears that aging is just a 

consequence of an extremely slow relaxation and is either interrupted for time scales larger 
than the (possibly astronomically large) equilibration time r eq or lasting forever within a 
spin glass phase with an infinite equilibration time. 

It has been suggested J7|,[| that aging is a manifestation of a slow domain growth at 
small temperatures, where after a certain waiting time t w a characteristic domain size R(t w ) 
is reached. The time dependence of subsequent observations (as for instance the thermore- 
manent magnetization) show a clear cross-over from dynamical processes characterized by 
length scales smaller than the already achieved domain size (and mainly taking place in- 
side the domains, for which reason it is called the quasi-equilibrium regime) to processes on 
larger time scales dominated by the continuation of domain growth through the movement 
of domain walls across the system. 

An activated dynamics scenario [0 is built upon the hypothesis that after a temperature- 
quench from a fully disordered state to low temperatures domains on a length scale L 
can grow by surmounting (free) energy barriers B(L) oc Lr, where ip < d — 1 is some 
(temperature-independent) exponent. Via thermal activation such a process needs a time 
t oc exp(AB(L)/T). Thus after a waiting time t w the characteristic domain-size is given by 

R(t w ) oc (logtj 1 ^ . (1) 

Originally such a scenario was proposed for the dynamics of a strongly disordered system 
well below its phase transition temperature. In the two-dimensional Ising spin glass, which 
does not have a finite temperature phase transition (see e.g. ||1U||), this scenario cannot hold 
for waiting times that are comparable to the (temperature-dependent) equilibration time r eq . 
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However, for t w <C r eq it might be valid at very low temperatures. Recent experiments on the 
two-dimensional short range Ising spin glass Ru2Cuo.89Coo.11F4 @,|] indicated the validity 
of an activated dynamics scenario. However, since the domain-size or the correlation length 
of spatial correlations has not been measured directly (e.g. via neutron scattering), the 
growth-law (|I]) had to be verified indirectly by establishing the activated scaling behavior 
I \oguj\/ \ogt w of the susceptibility Xtwi 1 ^) measured after aging the system for a waiting time 
t w . From this it was then concluded that ([!]) should also hold. 

On the other hand, in computer simulations one has immediate access to the quantities 
of interest and the validity of a logarithmic growth law like (H) can be checked directly. 
Therefore we present in this letter the results of extensive Monte-Carlo simulations of the 
aging scenario and the domain growth in the two-dimensional Ising spin glass model (more 
details will be published elsewhere ||11||). The latter is defined by the Hamiltonian 

H = - JijSiSj , (2) 

<ij> 

where Si = ±1, the sum extends over all nearest neighbor pairs of a L x L square lattice with 
periodic boundary conditions and interaction strengths J^- are Gaussian random variables 
with zero mean and variance one. We choose the usual heat bath algorithm for the spin 
dynamics and time is measure in Monte-Carlo sweeps through the whole lattice. We used 
lattice sizes L ranging from 30 to 200 depending on the temperature T (note that one does 
not need to consider system sizes much larger than the typical correlation length to exclude 
finite size effects). The quantities of interest are strongly fluctuating from sample to sample, 
for which reason we averaged over several thousand disorder realizations). The simulations 
were performed on a Parsytec GCell024 transputer cluster. 

A straightforward way to establish aging in the present model is to calculate the spin 
autocorrelation function 

tfftO = jj y E[(Si(t + t v ,)s i {t v ,))] m , (3) 

i 

where (• • •) means a thermal average (i.e. an average over different realizations of the thermal 
noise) and [• • -] av means an average over different realizations of the bond- disorder. For each 
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run the system is initialized in a random initial configuration corresponding to a quench 
from infinite temperature to the temperature T at which the simulation is done (e.g. S_(t w ) 
then denotes the configuration of the system at time t w after the quench). 

In figure 1 we show the results for C(t,t w ). For higher temperatures one observes that 
curves for waiting times t w larger than a particular value, which we denote with r eq , collapse. 
One says that aging is interrupted here |12|] and this is simply a manifestation of the fact that 



the system is equilibrated now and the spin correlations [{Si{t)Si{t'))\ aw become stationary 
in time (i.e. functions oit — t' only). For lower temperatures equilibration is not achieved on 
the time-scales explored and the resulting picture is indistinguishable from that observed in 
the three-dimensional Ising spin glass [13[] or the mean field Ising spin glass within the 



spin glass phase. A characteristic feature of the low temperature (T < 0.2) behavior is the 
crossover at t w from a very slow decay at times t < t w to a faster, algebraic decay at times 
t ^> t w . 

C(t,t w ) ocr« T) fort M ,<r eq , (4) 

which, however, is expected to fail at much larger time-scales t w 3> r eq . The exponent ((T) 
decreases with temperature and is e.g. 0.05 for T = 0.2. Although this is a rather small 
exponent our data exclude a logarithmic decay that is predicted for activated dynamics [|7| . 
In figure 2 we show that C(t,t w ) obeys simple dynamic scaling 

C(t,t w ) (xc[t/r(t w )] , (5) 

where the time scale r(t w ) is depicted in the insert. It is equal to t w for t w <C r eq and 
saturates (observable only for high enough temperatures) at r cq . Note that according to an 
activated dynamics scenario the scaling variable should be \og(t) / log(t w ) rather than t/t w 
(see e.g. equation 3.20 in 0). A t/t w scaling has also been observed in three-dimensional 
spin glasses fllBTJIB"! , the SK-model [[H]] and simplified spin glass models [|12|,|nj. 



The above described aging scenario can be thought of as being a consequence of a very 
slow domain growth, or equivalently, the slow increase of correlated volumes in the system 



after the temperature quench, as described in the introductory remarks. To verify this 
hypothesis we have to look at a quantity that provides us with the necessary informations. 
The domain growth in the various strongly disordered system like the site-diluted Ising 
model |I?fl , the random field Ising model |Tj| and the random bond ferromagnetic Ising 
model [TSj or the random bond Ising-chain EH] has already been investigated numerically. 



These models have the advantage that their ground state is known to be ferromagnetic, 
which makes the identification of domains easy. This is not the case for the present system 
- which is the reason why such an investigation is much more difficult here. One way to 
measure the domain size in spin glasses, which is however, very (computer-) time consuming, 
has been proposed in pJJ by studying the three-dimensional Ising spin glass. 
We decided instead to calculate an appropriate spatial correlation function 

G(r,t w ) = - £ [(S (t)S r (t))\ v , (6) 

Z w t =t w +l 

which is a generalization of the usual equilibrium correlation function G eq (r) = 
[{Trs SoS r exp — /3H(S.)} 2 } av to the present non-equilibrium situation (note that 
lirn^-^oo G(r, t w ) = G cq (r)). It is obvious that for small waiting times t w the correlations 
decay rapidly to zero on the length scale of one lattice spacing since the initial configuration 
is random and spins on different sites are uncorrelated. For increasing waiting times the 
system tries to relax into energetically more favorable configurations (closer to the — in this 
case unique but unknown — ground state). This is the process of domain growth mentioned 
in the beginning and can be read off from @ by an increase in the number of spins that are 
longer correlated (to contribute to the sum over t) over longer distances r. Thus the length 
scale of the decay of G(r,t w ) is a measure for the averaged domain size R(t w ). 

Usually when calculating G eq (r) the thermal average (■ • •) is replaced by a time average 
over a very long Monte-Carlo run. To get good statistics for the quantity (^j), especially 
for small waiting times t w , one has to average instead over many Monte-Carlo runs (the 
number of which decreases with increasing t w , fortunately) using different initial conditions 
and different thermal noise. Furthermore one performs a spatial average by taking into 



account all spin-pairs that are r sites apart in the x— as well as in the y-direction. To avoid 
a positive bias of the statistical errors induced by the square (which is necessary since the site 
correlations alternate in sign) we simulated two replicas a and b of the system simultaneously 
(both with identical realization of disorder, but with uncorrelated initial conditions and 
thermal noise) and calculated the four spin correlations (Sfi (t) (t) S^t) Sj.(t)) ■ This yields 
results identical to those obtained by using (So(t)S r (t)) 2 , as we checked explicitly. 

Figure 3 shows the result for G(r,t w ) for different waiting times t m at T = 0.8 One 
observes that the correlations grow steadily with increasing waiting time, until, at low enough 
temperature, the resulting curves collapse for t w > r eq . This is in full agreement with the 
picture that emerged from the time-dependent autocorrelation function C(t,t w ) described 
above. In a the linear-log plot of figure 3 the curves are only slightly bended, meaning that 
the decay is roughly exponential with a characteristic length scale. We define an effective 
correlation length (that is identical to this length scale for a pure exponential decay) via 

/•OO 

£(t w ) = 2 drG(r,t w ), (7) 
Jo 

which, as we mentioned above, is expected to be comparable to the averaged domain size 
R(t w ) P^ , P^| . The result is depicted in figure 4, where it can be seen that £(t w ) increases 
slowly on a logarithmic scale until it saturates at £ eq for t w > r eq . We did a least square fit 
to a logarithmic growth law 

£(t w ) -6 oc (\ogt w )W for^<r eq (8) 

and obtained ip — 0-63 ± 0.05 (roughly independent of temperature), which is in agreement 
with the bound ip < d — 1 and slightly smaller than the value found experimentally . This 
is compatible with an activated dynamics scenario, see eq. (|1]). However, one obtains also 
a good result (in terms of the outcome of a chi-square test) for a fit to an algebraic growth 
law 

£(t w )(xtZ iT) fort w <r cq . (9) 
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The exponent ct(T) decreases with temperature, being for instance 0.08 for T = 0.4 and 
0.04 for T = 0.2. These are rather small exponents, which is the reason why it is hard to 
discriminate between @ and (|9]). 

^From our point of view an algebraic domain growth yields a picture that is more con- 
sistent when taking into account our results for the autocorrelation function. On one side 
the proposition (H) implies a logarithmic scaling of the autocorrelation function C(t,t w ), 
which is not in agreement what we find, see (||), and also implies a logarithmic decay 
C(t,t w ) oc \og(t)~ y/ ^ for t > i w (< r cq ) 0, which disagrees with our result @. On the 
other hand these results (|J), (0) and @ can be explained consistently in the following way: 

Let us assume a modified activated dynamics scenario for domain growth in which 
the (free) energy barriers B(L) scale with the length scale L of the domains as B(L) ~ 
a(T) log(L) instead of B(L) oc and let us follow [|7| otherwise. Then one obtains for the 
domain-size R(t w ) oc ^ , which is exactly (§) with a(T) = a(T)A(T)/T. Furthermore 
from C(t,t w ) oc [R(t w ) / R(t)] x ' for t 3> t w (<^ r cq ) one concludes the scaling (|5|) as well as 
the algebraic time decay with ((T) = \'a(T). We should mention that an algebraic 
growth law like (^) for the domain sizes has been used by Koper and Hilhorst as a working 
hypothesis in their domain theory ||. 

In summary we have shown that aging is present at low temperatures even in the two- 
dimensional Ising spin glass model and is caused by a very slow domain growth. We have 
shown that the autocorrelation function scales like t/t w and decays algebraically for waiting 
times smaller than the equilibration time. Fitting the time dependence of the averaged 
correlation length or domain size to a logarithmic growth law or to an algebraic growth law 
yields reasonable results in both cases. However, an algebraic fit is slightly better and is 
also compatible with the results for C(t,t w ), whereas a logarithmic law is not. Results for 
the three-dimensional Ising spin glass |]13[ indicate that here the situation is quite similar. 



We think that the simultaneous investigation of correlations in time and in space are 
very important to obtain a consistent picture of the dynamical processes in spin glasses 
at low temperatures. It would be very helpful if experimentalists could find a practicable 



way to gain direct access to spatial correlation functions in spin glasses, as is, for instance, 
possible in random field systems ||24|| . In numerical simulations, however, such an endeavor 
is promising, as we have demonstrated here. For instance it would be of great interest to 
check such concepts as the overlap length in connection with temperature or field cycling 
experiments via a direct computation of the appropriate correlation functions in Monte- 
Carlo simulations (similar to those presented in this paper). That might help to resolve a 
long lasting debate among physicists [pt[25H27jl . 

We would like to thank the Center of Parallel Computing (ZPR) in Koln for the generous 
allocation of computing time on the transputer cluster Parsytec-GCell024. This work was 
performed within the SFB 341 Koln-Aachen-Julich. 
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FIGURES 

FIG. 1. Autocorrelation function C(t, t w ) as a function of time t for t w = 5 n (n = 1, . . . , 8) at 
T = 1.0 and 0.8, (n = 2, . . . , 8) at 0.6 and 0.2. The system size is L = 100 and the disorder average 
was performed over 256 samples. The errorbars are smaller than the symbols. 

FIG. 2. Scaling plot of the autocorrelation function C(t,t w ) versus t/r(t w ) at T = 0.2, The 
insert shows r(t w ). 

FIG. 3. The spatial correlation function G(r,t w ) for t w = 5 n (n = 1, . . . , 8) at temperatures 
T = 0.8. The insert shows a scaling plot of G(t,t w ) versus r/p(t w ). The scaling length | is equal 
(within the errorbars) to £(t w ) defined in the text and shown in the next figure. The system size 
is L = 100 and the disorder average was performed over 256 samples. 

FIG. 4. Log-linear plot of the correlation length £(t w ) as defined in the text. The full lines are 
least square fits to (log^) 1 /^ with V = 0.66 for T = 0.2, V = 0.60 for T = 0.3 and tp = 0.63 for 
T = 0.4. The insert shows £(t w ) at somewhat higher temperatures. 
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